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We present the first results for anisotropy searches around the galactic center at EeV energies using data from 
the Pierre Auger Observatory. Our analysis, based on a substantially larger data set, do not support previous 
claim of anisotropy found in this region by the AGASA and Sugar experiment. Furthermore we place un upper 
limit on a possible point like source located at the galactic center which exclude several scenarios predicting 
neutron sources in this location. 



1. Introduction 

The Galactic Centre (GC)is an attractive tar- 
get for the studies of cosmic ray (CR) anisotropy 
at EeV energies. The GC contains objects that 
might be candidates for powerful CR accelerators, 
and the recent high significance observation by 
H.E.S.S. of a TeV 7 ray source near the location 
of Sagittarius A* pQ , together with the discovery 
of extended emission from giant molecular clouds 
in the central 200 pc of the Milky Way [2], mo- 
tivates the search for excesses in this direction. 
The GC passes only 6° from the zenith at the 
site of the Pierre Auger Observatory in the south- 
ern hemisphere and makes it particularly suitable 
for anisotropy studies in this region. The num- 
ber of CRs of EeV energies accumulated so far 
at the Pierre Auger Observatory from this part 
of the sky greatly exceeds that from previous ob- 
servations. We adopt hereafter for the GC co- 
ordinates the Sagittarius A* J2000.0 coordinates, 
(a, 5) = (266.3°, -29.0°), 

There have also been reports by the AGASA 
experiment lj indicating a 4.5er excess of cos- 
mic rays in the energy range 10 18 -10 18 4 eV and 
by the Sugar experiment |Hj • Note that the GC 
itself lies outside the AGASA field of view (6 > 
-24.2°). 

We also searched for a point-like excess from 
the GC, as EeV neutrons emitted by one of the 
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possible energetic sources in the centre of the 
Galaxy may reach the Earth before decaying, 

After a brief description of the Auger detector 
we discuss the angular resolution of our surface 
array and the sky coverage estimation. We fol- 
low with a presentation of the data set used in 
this analysis and with the results concerning the 
AGASA and SUGAR claims as well as an upper 
limit on a neutron source located at the centre of 
our galaxy. 

2. The Auger detector. The hybrid design 

A fundamental characteristics of the Pierre 
Auger Observatory{5] is its capability of hybrid 
reconstruction of cosmic ray showers 7 . Two in- 
dependent detectors, the Surface Detector (SD), 
which samples the shower particles at ground, 
and the Fluorescence Detector (FD), which col- 
lects the fluorescence light emitted by the shower 
particles along their path in the atmosphere, are 
able to measure the energy and direction of the 
same cosmic ray shower. The enhanced capabili- 
ties of the Auger hybrid detector are examined in 
this paper. 

The Pierre Auger Observatory was designed to 
observe, in coincidence, the shower particles at 
ground and the associated fluorescence light gen- 
erated in the atmosphere. This is achieved with 
a large array of water Cherenkov detectors cou- 
pled with air-fluorescence detectors that overlook 
the surface array. It is not simply a dual exper- 
iment. Apart from important cross-checks and 
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measurement redundancy, the two techniques see 
air showers in complementary ways. A single air 
shower is detected 3-dimensionally. The ground 
array measures the 2-dimensional lateral struc- 
ture of the shower at ground level, with some 
ability to separate the electromagnetic and muon 
components. The fluorescence detector records 
the longitudinal profile of the shower during its 
development through the atmosphere. 

A hybrid event is an air shower that is simul- 
taneously detected by the fluorescence detector 
and the ground array. The Observatory was orig- 
inally designed and is currently being built with 
a cross-triggering capability. Data are recovered 
from both detectors whenever either system is 
triggered. If an air shower independently triggers 
both detectors the event is tagged accordingly. 
There are also cases where the fluorescence detec- 
tor, having a lower energy threshold, promotes a 
sub-threshold array trigger. Surface stations are 
then matched by timing and location. This is an 
important capability because these sub-threshold 
hybrid events would not have triggered the array 
otherwise. The geometrical reconstruction of the 
air shower's axis is accomplished by minimizing a 
X 2 function involving data from all triggered ele- 
ments in the eye and on the ground. The recon- 
struction accuracy is far better than the ground 
array counters and the single eye could achieve 
independently 0. 

The combination of the air fluorescence mea- 
surement and the particle detection on the ground 
provides an absolute energy calibration. 

3. Surface Array Angular Resolution 

The angular resolution for the SD is deter- 
mined, on an event by event basis, from the zenith 
(9) and azimuth (0) uncertainties obtained from 
the geometrical reconstruction, using the relation: 

F(r)) = 1/2 [V[9] + sm 2 (9) V[<f>]) (1) 

where n is the space- angle, and V[9] and V[4>] are 
the variance of 6 and <j> respectively. If 9 and 
<fi/sxn(9) have Gaussian distribution with vari- 
ance <7 2 , then F(rj) = a 1 and r\ has a distribution 
proportional to e~ v l 2a d(cos(rj))d<j). Then, if we 
define the angular resolution (AR) as the angular 



radius that would contain 68% of showers coming 
from a point source, AR = 1.5 \/F(rj). 

The angular resolution depends strongly on the 
timing resolution of the water Cherenkov detec- 
tors (WCDs) and weakly on the shower front 
model and the core position uncertainty. The 
WCDs timing uncertainty is directly modeled 
from the data fsection l3.1[l . This model is based 
on the physics of the shower and the measurement 
process. It can be adjusted using two pairs of ad- 
jacent stations located in the surface array (sec- 
tion 13. 2|) . The model is validated by studying the 
X 2 probability distribution for the geometrical re- 
construction (section !^. 3|l . The angular resolution 
is estimated for the SD-only reconstruction and 
by comparison with the hybrid data fsection l3.4fl . 

3.1. The Time Variance Model 

The angular accuracy of the SD events is driven 
by the accuracy with which one can measure the 
arrival time (T s ) of the shower front in each sta- 
tion. The particle arrival time in the shower front 
can be described as a Poisson process over some 
interval time T. The first particle arrival time is 
used as the estimator for the shower front arrival. 
It is given by 2 T\ = T s +ti, where T s is the shower 
front time and t\ is taken to follow an exponential 
distribution function with decay parameter r. 

Since we estimate the parameter T from the 
data itself, the variance of T± (given by the vari- 
ance of t\) is: 

The variance of the arrival of the first particle 
in the SD stations, taking into consideration the 
GPS uncertainty and the resolution of the flash 
analog-to-digital converters (FADC), can then be 
written as: 

V[T 1 ] = a* (£=J»Y n —±+b 2 (3) 
\ n J n+1 

where T50 is the time interval that contains the 
first 50% of the total signal as measured by the 

2 In fact, an unbiased estimator should be To = Ti — E[ti], 
where E[ti] is the expectation value. 
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photomultiplier FADC traces. The two free pa- 
rameters a and b can be determined with the adja- 
cent station data. We expect that the parameter 
a should be close to 1, while b should be given 
by the GPS clock accuracy (about 10 ns) and the 
FADC trace resolution 25/VT2 ns, that is b ~ 12 
ns. 

To calculate the number of particles (n) we as- 
sume that all particles hit the detector with the 
same direction than the shower axis, and that the 
muons are mostly the ones that contribute to the 
time measurements. Then, we obtain n as the ra- 
tio between the total signal (S) in the WCD and 
the average track length, TL(6), of the particles. 

The average track length can be computed as 
the ratio of the detector volume (V) and the area 
(A) subtended by the arriving particles, and is: 



V 

TL{9) = - 



7rr 2 cos(0) + 2rh$m(6) 



(4) 



where is the zenith angle, r = 1.8 m is the detec- 
tor radius, and h = 1.2 m is the detector height. 

3.2. Testing the model with doublets 

Two pairs of adjacent surface detector stations 
separated by 11 m ("doublets") have been in- 
stalled in the field of the Auger Observatory. 
These pairs enable comparison of timing and sig- 
nal accuracy measurements. We used the data 
of the doublets to verify the time variance model 
and also to adjust the constants a and b from it. 
For each event we computed the time difference 
as AT = dT« - dT& where <2T« (cff( 2 )) is the 
time difference from the first (second) detector of 
the doublet to the fitted shower front. Doing that, 
AT does not depend on the shower front shape, 
since the twin detectors are very close to each 
other. We used 1693 events (from April/2004 to 
June/2006) to fit for the two parameters a and b, 
and we obtained: 

a 2 = 0.98 ± 0.05 
b 2 = 150 ns 2 ± 18 ns 2 

which is in good agreement with our expectations 
(a 2 = 1, b 2 = 144 ns 2 ). 
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Figure 1. T he RMS of the distribution of 
AT/ y/V[AT], as a function of the shower zenith 
angle (top), the average signal in the doublet de- 
tectors (middle), and the distance to the shower 
core (bottom). 



3.3. Validation of the Time Variance 
Model 

If the time variance model describes 
correctly the measurement uncertainties, 
the distribution of AT/y/V[AT], where 



V[AT] 



V[T^} 



V[T^}, should have 
unit variance. In figure ^ we show the RMS of 
the distribution of AT/ y/V[AT] for the doublets 
as a function of cos((9) (top), the average signal 
(middle), and the distance to the core position 
(bottom). In all the cases, the RMS is almost 
constant and close to unity. 

In figure we plot the x 2 probability distribu- 
tion from our geometrical reconstruction proce- 
dure using as timing errors the model above. We 
only plot probabilities larger than 1% to avoid 
the large peak at zero corresponding to badly re- 
constructed events (~ 9%). These distribution 
are almost flat indicating that our error model is 
correct. 

3.4. Surface Detector Only 

Given the correctness of our error model, we 
can extract the angular resolution on an event 
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Figure 2. The x 2 probability distribution for 
all events (top), events with zenith angle smaller 
than 55° (middle), and events with zenith angle 
larger than 55° (bottom). In the last figure the 
distribution is plotted with two different scales, 
the same than the others (full line) for compari- 
son reasons and a zoom (dashed line) to see the 
details. 



by event basis directly out of the minimization 
procedure. 

The angular resolution is about 2.7° in the 
worst case of vertical showers with only 3 stations 
hit. This value improves significantly for 4 and 5 
stations 3 . For 6 or more stations, which corre- 
sponds to events with energies above 10 EeV, the 
angular resolution is in all cases better than 1°. 

All quoted errors are statistical only. We did 
not, at this stage, investigate possible biases or 
systematics in the determination of the arrival di- 
rection angles. 

3.5. Comparison with Hybrid events 

Finally, in figure 01 we show the space angle 
between the SD-only and hybrid geometrical re- 
constructions for showers with different number 



2.5 - 



& 2 - 



+ 



1.5 



3 4 and 5 stations events have the same number of degrees 
of freedom in the geometrical fit, hence similar resolution. 
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Figure 3. Angular resolution (AR) for the SD as 
a function of the zenith angle (9) extracted from 1 
reconstruction procedure. The AR is plotted for 
various stations multiplicities (circles: 3 stations, 
squares: 4 stations, up triangles: 5 stations, and 
down triangles: 6 or more stations). 



of stations and different zenith angle ranges. The 
distributions plotted were fitted with a Gaussian 
resolution function dp oc e~ n / 2<T d(cos(r]))d<j), 
where r\ is the space-angle. The a obtained in 
the fit is related to the angular resolution by 
RA = 1.5 a. 

The parameter r\ given by the fit is in good 
agreement with the value from figurc|3| using an 
hybrid resolution parameter r] of 0.5 — 0.6°. 

The angular resolution of the surface detec- 
tor is then found to be better than 2.7° for 3- 
fold events, better than 1.7° for 4-fold and 5- 
fold events and better than 1.0° for higher multi- 
plicity (which corresponds to energies larger than 
10 EeV). 

4. Coverage 

To study anisotropy, one needs the background 
expectations for the different sky directions under 
the assumption of an isotropic CR distribution. 
Modulations of the exposure in right ascension 
are induced by the dead time of the detectors and 
by the constantly growing array size. So those ex- 
perimental issues (that are carefully recorded by 
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Figure 4. Comparison between hybrid and SD- 
only geometrical reconstructions. Top, for 3 sta- 
tions with two zenith angle ranges (0° < < 30° 
and 30° < 6 < 50°). Middle 4 stations (left) 
and 5 stations (right) with 30° < 6 < 50°. Bot- 
tom, for 6 stations or more with two zenith angle 
ranges (30° < 6 < 50° and 50° < 6 < 50°). 



the auger data acquisition system) must be fully 
taken into account. Also weather variations, may 
be non-negligible since they affect the shower de- 
velopment in the atmosphere and/or the response 
of the electronics. 

Preliminary studies of these effects have shown 
that the possible weather-induced background 
modulations for the present data set are at a level 
of 1%. This is below the Poisson noise for the an- 
gular windows considered 4 . 

We have followed two different technics JT] to 
estimate the coverage for the SD analysis: 

• The semi-analytic technique: The 

zenith angle distribution of the events in the 
considered energy range is fitted and then 
convoluted with the number of hexagons 



4 A detailed account of weather effects is certainly neces- 
sary to test large scale patterns at the few percent level. 
Relevant studies are in progress. 



with active detectors (which gives a mea- 
sure of the aperture for events satisfying 
the quality trigger criterion) as a function 
of time, assuming a uniform response in 
azimuth. Through this procedure one ob- 
tains an exposure which accounts for the 
non-saturated acceptance effects and for 
the non-uniform running times and array 
growth. 

• The shuffling technique: Here the ex- 
pected number of events in any direction 
is obtained by averaging many data sets 
obtained by shuffling the observed events 
in the energy range of interest so that the 
arrival times are exchanged among them 
and the azimuths are drawn uniformly. 
The shuffling can be performed in separate 
zenith angle bins or by just mixing them 
all, and we found no significant difference 
between these two possibilities. By con- 
struction, this exposure preserves exactly 
the 9 distribution of the events and ac- 
counts for the detector dead times, array 
growth and even for weather-induced mod- 
ulations. It might however partially absorb 
modulations induced by large scale intrinsic 
anisotropy present in the CR flux, such as 
those due to a global dipole. 

The background estimate obtained with the 
two technique differ by less than 1% in the GC 
region. It is 0.5% larger with shuffling than the 
semi-analytic method. This difference is much 
smaller than the size of the excesses we are test- 
ing and also below the Poisson fluctuations, In 
the following we mainly quote the values obtained 
using the semi-analytic technique. 

5. Data set 

The Auger surface detector [Jj], has been grow- 
ing in size during the data taking period consid- 
ered in this work. 154 detectors deployed in Jan- 
uary 2004 and 930 in March 30 th 2006 

The basic cell of surface detectors is triangu- 
lar, with separations of 1.5 km between detector 
units, and hence the complete array with 1600 
detectors will cover an area of 3000 km 2 . 
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We consider the events from the surface detec- 
tor (SD) array with three or more tanks triggered 
in a compact configuration. The events have to 
satisfy the level 5 quality trigger condition, which 
requires that the detector with the highest signal 
be surrounded by a hexagon of working detec- 
tors, since this ensures that the event is well re- 
constructed. We also restrict the events to zenith 
angles 9 < 60°. 

The energies are obtained using the inferred 
signal size at 1000 m from the reconstructed 
shower core, 5(1000), adopting a conversion that 
leads to a constant flux in different sky directions 
above 3 EeV, where the acceptance is saturated. 
This is the so-called Constant Intensity Cut crite- 
rion implemented in |12j . A calibration of the en- 
ergies is performed using clean fluorescence data. 
The estimated systematic uncertainty in the re- 
constructed shower energy with the fluorescence 
technique is currently 25% [T5) . 

From the uncertainty in the measurements of 
the signals from the Cherenkov tanks JSj the sta- 
tistical uncertainty in the energy determination 
which results from the fitting procedure is about 
20% for the energy range considered in this work, 
i.e. 10 17 - 9 eV < E < 10 18 5 eV. Notice that in 
this energy range 48% of the events involve just 
three tanks, 34% involve 4 tanks and only 18% 
more than 4 tanks. For three tank events the 
68% quantile angular resolution is about 2.2° and 
the resolution improves for events with 4 tanks or 
more 

Regarding the hybrid events, i.e. those with 
signal from both the fluorescence detectors (FD) 
and surface array, the angular resolution achieved 
is much smaller, typically below 1 degree 
Also, given that hybrid events may trigger with 
just one surface detector, the associated energy 
threshold (~ 10 17 eV) is lower, and events up to 
zenith angles of 75° are included in the data set. 
However, the statistics accumulated are signifi- 
cantly less, in part due to the ~ 15% duty cycle 
of the fluorescence telescopes and also because at 
EeV energies the FD is not fully efficient at de- 
tecting showers over the full SD array. 

There are for instance 79265 SD events in the 
data set considered with energies 10 17 9 eV < 
E < 10 185 eV, while the corresponding number 



of well reconstructed hybrid events in the same 
energy range is just 3439. Note that ~ 25% of 
the hybrid events in this energy range involve less 
than three surface detectors, and are hence not 
included in the SD only data set. 

6. Results 

All results presented in this section are part of 
a dedicated publication^] 

7. AGASA and SUGAR excesses 

The map of overdensities around the GC re- 
gion is shown on Figure [S] Significances in this 
map are calculated in circular windows of 5.5° us- 
ing a Gaussian approximation to the Poisson law 
(counts in the 5.5° windows are well above 100) 
The angular scale chosen and the energy range 
shown (10 17 ' 9 -10 18 5 eV) mimics the Sugar search 
window and is also convenient to visualize the dis- 
tribution of overdensities in the window explored 
by AGASA. The galactic plane is shown as a solid 
line (a cross for the GC) The big circle region in 
which AGASA reported an excess (in a slightly 
narrower energy range). The dashed line indi- 
cates the lower boundary of their observations. 

Overdensities present in this map are consistent 
with the expectation from statistical fluctuations 
of an isotropic sky. On Figure we show the 
distribution of these overdensities together with 
the expectations from an isotropic flux and no 
significant departure from isotropy is observed. 

In the 20° circle centered at the AGASA loca- 
tion and for energies between 10 18 eV < E < 
10 18 4 eV, we counted 2116 events while while 
2169.7 are expected using the shuffling technique 
( 2159.6 using the semi-analytic) No significant 
excess is observed. 

There may be systematic differences in the en- 
ergy calibration of the two experiments. To test if 
these differences could possibly mask the AGASA 
excess, we calculated the observed and expected 
rates for different energy ranges in offsets of 0.1 
decade keeping E max /E min fixed. We have also 
added a systematic error of 1% to the expected 
rates to account for possible weather induced 
modulations. We found no significant excesses 



Figure 5. Map of CR overdensity significances near the GC region on top-hat windows of 5.5° radius. 
The GC location is indicated with a cross, the galactic plane as solid line. 



in the AG AS A region for any of these cases. In 
particular, at the 2a level (95% CL) the excess in 
this region is always less than 6%, well below the 
22% excess reported by AGASA. 

The acceptance of Auger in this energy range 
is not yet saturated, it is larger for heavy nu- 
clei than for protons because showers initiated by 
heavier primaries develop earlier and are hence 
more spread out at ground level. Using ref [T9"| 
for the acceptance of p and Fe primaries, we find 
that the sensitivity to protons is about ~ 30% 
smaller than to Fe in the energy range studied 
(assuming an E~ 3 spectrum). If the 22% excess 
reported by AGASA (which had full efficiency at 
EeV energies) was due to nucleons while the back- 



ground was due to heavy nuclei, at least a 15% 
excess should have been expected in Auger data. 
This is much larger than the upper limit we are 
obtaining. 

For the excess reported by the SUGAR ex- 
periment, we find in the same angular window 
and energy range that n b s /n exp — 286/289.7 = 
0.98 ±0.06, and hence with 10 times the statistics 
we found no no significant excess in this window. 

8. Neutron source at the GC 

8.1. Surface detector 

The energy range for this search is between 
E m in = 10 17 - 9 eV and E max = 10 18 5 eV. Be- 
low E m i n the Auger SD acceptance is suppressed 
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Figure 6. Histogram of overdensities on 5.5° ra- 
dius windows and for 10 17 9 eV < E < 10 18 5 eV, 
together with the average isotropic expectations 
(gaussian of unit width and zero mean). The 
histrogram is build over the same sky patch as fig- 
ure 6 (-65° < dec < 5° and 160° < RA < 360°). 



and most neutrons from a source at the GC would 
have decayed before reaching the Earth. Energies 
above E max are assumed to be be hard to achieve 
for galactic sources. 

The optimal search window for a point-like 
source is obtained using a Gaussian filter match- 
ing the angular resolution of the experiment [20] . 
Following section we used a 68% quantile a of 
2.25° for the average resolution of the events in 
this energy range (in fact more than 80% of the 
events have a better resolution than this). 

In the direction of Sagittarius A* we get 
n obs /n exp = 53.8/45.8 (a ratio of 1.17 ± 0.10). 
Applying the results of [201, we get a 95% CL 
upper bound on the number of events from the 
source of n 9 s 5 — 18.5. 

Assuming : 

• that the source spectrum shape is similar to 
that of the bulk CRs, (cx E~ 3 in this energy 
range), 



• that both the bulk CRs and the source CR 
are protons in this energy range (We discuss 
how the limit is modified if the bulk CRs 
were heavier) , 



we have [20] : 



(5) 



Where §cr and § 9 S 5 are the bulk CR flux and 
the source 95% flux upper limit integrated in the 
energy range under study and where n exp and n 9 , 5 
are as above. 77 is the parameter of the optimal 
Gaussian window that corresponds to our angular 
resolution 5 (here r\ = 2.25/1.5 = 1.5°) 

For the differential CR spectrum flux we take 



§cr(E) ~ i 30 



E 
EeV 



EeV^km-V^sr -1 : 



(6) 



which has an E~ 3 dependence. The factor £ is 
close to unity and parametrizes the uncertainties 
in the CR flux normalization, so that the flux 
bounds will be simply proportional to £. 

From eq. El we get a 95% CL upper bound for 
the source flux integrated over the energy range 
considered of : 



§f = £ 0.08 km"^ yr" 



(7) 



If the bulk of the CRs were heavy nuclei the 
upper limit to the flux from the source need to to 
be scaled by a factor ~ 1.3 under the assumption 
that the CRs are iron nuclei and that the source 
is of neutrons. We thus see that the bound on the 
neutron flux could be up to ~ 30% higher if the 
CR composition at EeV energies were heavy. 

Theoretical predictions for neutron fluxes (only 
those associated with the AGASA claim) are 
based on the AGASA normalization. This nor- 
malization is about a factor of 3 larger than 
the Auger one so these predictions must be re- 
duced by this factor to be compared with our 
flux bound. The predictions of refs. [23, [23 and 
|23j . which exceed the upper-bound obtained by 
more than one order of magnitude, are already 
excluded, and that of [23 is at the level of the 
present Auger sensitivity. 

5 Remember that rj and the 68% resolution quantile a are 
related by a = l.5rj (cf section 151. 
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8.2. Hybrids 

As discussed earlier, hybrid events, detected by 
both the FD and SD have a better angular resolu- 
tion, 0.7° at 68% C.L. in the energy range studied. 

For energy between 10 17 9 eV < E < 10 18 5 eV, 
no significant excess is seen in the GC direction. 
Using, for instance an optimal top-hat window 
of 1.59ct ~ 0.75° radius, 0.3 events are expected 
(shuffling technique estimate) while no event di- 
rection falls within that circle. 

This gives to a source flux upper-bound at 95% 
CL of 

$f = £0.15knr 2 yr-\ (8) 

which is about a factor of two weaker than the 
SD flux bound. 

Note that quality requirement for hybrid events 
is to have the maximum of the shower develop- 
ment inside the field of view of the telescopes, this 
affects the sensitivity to different primaries. The 
bound obtained here is a conservative one if the 
bulk of the CRs are heavy nuclei. 

9. Conclusions 

The Pierre Auger Observatory is a hybrid de- 
tector with excellent capabilities for studying the 
highest energy cosmic rays. Much of its capability 
stems from the accurate geometric reconstruction 
it achieves. 

The construction of the Southern Observatory 
is well under way. Eighteen FD telescopes and 
more than 60% of the surface array are in oper- 
ation taking data routinely. At the present rate 
of deployment, construction will be finish in mid 
2007. Detectors are performing very well and the 
first results are very encouraging. 

The angular resolution of the surface detector 
was found to be better than 2.7° for 3-fold events, 
better than 1.7° for 4- fold and 5-fold events and 
better than f.0° for higher multiplicity (which 
corresponds to energies larger than 10 EeV). 

Using a bit more than the 2 first year of Auger 
data we have searched for localized anisotropy 
near the direction of the Galactic Centre. With 
statistics much greater than those of previous ex- 
periments, we have looked for a point-like source 
in the direction of Sagittarius A, without finding 



a significant excess. 

We exclude several scenarios of neutron sources 
in the GC as suggested recently. Our searches 
do not support the large excesses reported in 
AGASA data (of 22% on 20° scales) and SUGAR 
data (of 85% on 5.5° scales). 
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